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Abstract 

We propose that the fermionic superpartner of a weak-scale Goldstone boson can be a nat- 
ural WIMP candidate. The p-wave annihilation of this 'Goldstone fermion' into pairs of Gold- 
stone bosons automatically generates the correct relic abundance, whereas the XENONIOO 
direct detection bounds are evaded due to suppressed couplings to the Standard Model. 
Further, it is able to avoid indirect detection constraints because the relevant s-wave annihi- 
lations are small. The interactions of the Goldstone supermultiplet can induce non-standard 
Higgs decays and novel collider phenomenology. 



1 Introduction 

Cosmological observations now provide overwhelming evidence that about 20% of the energy 
density of the universe is some unknown form of cold dark matter [1| . The most popular candidates 
are weakly interacting massive particles (WIMPs) which can produce the correct relic abundance 
after freeze out, 

^DM/i'-O.l/^. (1.1) 

{av) 

A natural candidate for WIMP dark matter arises in extensions of the Standard Model with low- 
scale supersymmetry (SUSY) and i?-parity. In such models the lightest supersymmetric particle 
(LSP) is automatically stable and generically has mass on the order of the weak scale [2]. 

The 'WIMP miracle' is the statement that a particle with a mass and annihilation cross section 
typical of the weak scale will automatically yield a relic abundance that is within a few orders of 
magnitude of the observed value. This paradigm has been challenged by recent direct detection 
searches for WIMPs. In particular, XENONIOO recently set the most stringent upper limit on 
the spin-independent elastic WIMP-nucleon scattering cross section, asi = 7.0 x 10~^^ cm^ = 
7.0 X 10~^ pb, for a 50 GeV WIMP at 90% confidence [3]. This large discrepancy between the 
necessary annihilation cross section and the direct detection bound is increasingly difficult to 
explain in the usual WIMP dark matter scenarios. 

For example, within the minimal supersymmetric Standard Model (MSSM), one must typi- 
cally tune parameters in order to explain this difference [1]. A standard approach is to consider 



parameters in which asi is suppressed below direct detection constraints. At generic points in 
the parameter space this will also imply a suppressed annihilation cross section and thus a relic 
abundance that is too large. In order to overcome this problem one needs to assume special rela- 
tions among a priori unrelated parameters in order to boost the annihilation rate. For example, a 
pure bino LSP would require coannihilation (due to an accidental slepton degeneracy) or resonant 
annihilation to obtain the correct annihilation cross section [5J. Alternately, the observation that 
Higgsinos and winos have annihilation cross sections that are typically too large allows one to 
tune the LSP to be a specific combination of bino, Higgsino, and wino to generate the correct 
abundance [6]. This 'well-tempered neutralino' scenario, however, is now strongly disfavored by 
XENONIOO H]. 

In light of this tension, it is natural to consider non-minimal SUSY models in which 

• the WIMP is a weak scale LSP, 

• the direct detection cross section is suppressed while maintaining the correct relic abundance 
without any fine tuning, and 

• the experimental prospects in near future include novel collider signatures. 

We therefore extend the MSSM by a new sector with an approximate global symmetry which is 
spontaneously broken in the supersymmetric limit. A natural WIMP candidate that satisfies the 
above criteria is the fermionic partner of the Goldstone boson which we refer to as the Goldstone 
fermion, x- This particle can naturally sit at the bottom of the spectrum because it lives in the 
same chiral supermultiplet as the Goldstone boson a and is thus protected by Goldstone's theorem 
and SUSY. Even when SUSY is broken, the Goldstone fermion can remain light with mass at or 
below MgusY [3 El E]- This scenario is a weak scale version of axino dark matter [TU]; for an 
early attempt containing similar elements see [TTj. Similar realizations also appear in dark matter 
models where the LSP has a large "singlino" component [12]; such models can reproduce the mass 
spectrum of Goldstone fermion dark matter but do not have a limit where the global symmetry is 
broken while SUSY is exact. In particular the singlino dark matter effective interactions do not 



come from an effective low-energy Kahler potential as discussed in Section 2.1 Further, due to 
singlino-Higgsino mixing, such models typically require tuning to avoid direct detection bounds. 

The SUSY non-linear sigma model is a generic low-energy theory of the Goldstone supermul- 
tiplet based only on the symmetry breaking pattern [13]. It can be organized as an expansion in 
inverse powers of the symmetry breaking scale, /. In particular the leading order contribution 
to dark matter annihilation is controlled by a trilinear derivative coupling Xl^l^xdi^a/ f ■ If the 
global symmetry is anomalous with respect the SM gauge group, the Goldstone bosons will, in 
turn, decay to stable SM particles, a ^ gg , 77. All the interactions can be perturbative and com- 
patible with gauge coupling unification if the mediators of the anomaly come in complete GUT 
multiplets. If the Goldstone fermion mass is around the weak scale and the symmetry breaking 
scale / is around the TeV scale, then the resulting annihilation cross section is automatically in 
the thermal WIMP range 

M^K/r)(T^/mJ^lpb. (1.2) 

The freeze-out temperature T/-/m^ ~ 1/20 is insensitive to details of the model and appears 
because XX aa is a p-wave process. 
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After electroweak symmetry breaking at fEw = 175 GeV, the CP-even scalar component of 
the Goldstone chiral multiplet mixes with the Higgs boson and generates an effective hxx couphng 
which is suppressed by m^VEw / f'^ ~ 0.01. While standard Higgsino-like dark matter in the MSSM 
gives a large direct detection cross section, Goldstone fermion scattering off nuclei lies just below 
the XENONIOO bound, 

asi ^ ' ^T"" ^ 10-^^ cm^ . (1.3) 

Note that this suppression factor in asi is roughly of the same order as the suppression needed in 
the annihilation cross section for a standard weak-scale WIMP, {<7v)y/iMP ~ ^Q^weak/(100 GeV)^ 
~ 150 pb, to obtain the correct abundance (1.1). 

Finally, Goldstone fermion dark matter has novel consequences on Higgs phenomenology at 
the LHC. The global symmetry requires a derivative coupling between the Goldstone boson and 
the Higgs boson ~ t>Ew//^(<9a)^/i. If kinematically allowed, the Higgs boson decays into four light 
unflavored jets, /i — )■ 2a — >■ 4j, with a sizeable branching ratio. This decay mode is 'buried' under 
the QCD background. Such non-standard Higgs decays have recently been investigated in SUSY 
models motivated by the little hierarchy problem [HI [TS]. For Goldstone fermion dark matter, 
the Higgs might only be 'partially buried' with a branching ratio of 30% to the Standard Model. 
Alternately, one can hope to discover the Goldstone boson itself by looking for the a ^ 2g decay. 
Together with the direct detection of its fermionic superpartner, such a discovery would be strong 
evidence that the dark matter particle emerges because of the Goldstone mechanism and SUSY. 

The paper is organized as follows. We introduce the effective low-energy theory of a Goldstone 
supermultiplet in Section [2] and extend this by including SUSY and explicit global symmetry 
breaking in Sections |3] and |4j Readers who are primarily interested in dark matter phenomenol- 
ogy can proceed directly to Sections |5] and [6| where we review (in) direct detection prospects and 
calculate the relic abundance. We discuss the LHC phenomenology in Section [7j In Appendix |A] 
we present simple models that realize this scenario. Details of the annihilation cross section cal- 
culation are given in Appendix [Bj Remarks on a possible Sommerfeld enhancement are presented 
in Appendix [C] 



2 The Goldstone Supermultiplet 

We consider a supersymmetric gauge theory with a global U(l) symmetry that is broken by fields 
"^i which obtain vacuum expectation values (vevs) /«. In the limit of unbroken supersymmetry, 
the theory has a massless Goldstone chiral superfield, 

A= ^{s + ia) + V2ex + e^F. (2.1) 
v2 

which is the low-energy degree of freedom of the high-energy fields, 

^^ = /^e'^'^/^ (2.2) 
where the effective symmetry breaking scale is 

P = Y.^fh (2.3) 
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and Qi is the U(l) charge of We refer to the component fields as the Goldstone boson a, 

the sGoldstone s, and the Goldstone fermion x- In models where the U(l) is a Peccei-Quinn 
symmetry, these are typically called the axion, saxion, and axino, respectively. The mass of the 
CP-odd scalar a is directly protected by the Goldstone theorem while the s and d masses are, in 
turn, protected by supersymmetry. 

The Goldstone boson shift symmetry acts on the chiral superfield as A ^ A + icf. It is thus 
often convenient to consider a non-linear realization of the Goldstone chiral superfield, G = e^^^ , 
which naturally transforms under the U(l) shift symmetry, G — > e^'^G. In the absence of explicit 
global symmetry breaking, this shift symmetry forbids any superpotential term involving A. 



2.1 Effective Kahler potential 

The shift symmetry restricts the dependence of the Kahler potential on the Goldstone superfield 
to take the form 

K ^ K{A + A\^j^). (2.4) 

We have written $l to denote light fields which are uncharged under the global symmetry. Note 
this general form includes the canonical term AA^ which is {A+A^y up to a Kahler transformation. 

We may examine the Goldstone self-interactions by expanding the canonically normalized 
Kahler metric in inverse powers of the scale /: 

^(2) = =l + hj{A + A^) + h^^{A + A^r + . . . , (2.5) 

where q is an reference U(l) charge of the theory. The choice of q is arbitrary and irrelevant since 
the combination f/q is invariant under charge rescaling. For simplicity we set g = 1 henceforth. 
After integrating out the auxiliary fields, the general form of the Lagrangian is 

- ^i^(3)(«) ix^^'xd.a) + I (k^,){s) - j (XX)(X^X^) , (2-6) 

where K(^n) — d'^K/dA^. Passing to four-component Dirac spinors and expanding the Lagrangian 
in inverse powers of 1// yields, 

1 + + 62-^^2 + ... j (^-d^^sd^s + \d^ad^a + '-xi^d^x^ (2.7) 

^ir2^1^''T'^") ^^^'^'^^ + (^2 - 6? + • • • ) [(XX)^ - [Xl'xf] . 

The coefficients 61^2, •■ completely characterize the self-interactions of the Goldstone multiplet in 
the symmetric limit. The hi coefficient is particularly important for the dark matter abundance 
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since it controls the size of the XX(^ vertex. The tree-level contribution to bi can be determined 
by comparing (2.5) to the canonical Kahler potential of the high-energy fields ipi, 



K = J2 ^1^^ = E /^^'(^+^')/^. (2.8) 

i i 

Note that in the absence of explicit U(l)-breaking terms, A does not get a vev and K is canonically 
normalized with respect to the Goldstone superfield. All Goldstone self-interactions are calculable 
from the physical Kahler metric, 

ka^' =j,Y: /^^-(^+^^)/^ = 1 + ^-^±^ (^y: ^'f^) +■■■■ (2-9) 

In particular, the tree-level contribution to the bi coefficient is given by 

^i = 4E^'/'- (2.10) 

Note that bi is invariant under overall charge scaling. In simple models with just two fields of 
opposite charge, bi is bounded, — 1 < &i < 1. In general, however, there is no such restriction on 
bi in theories with more fields or with dynamical U(l) breaking. 



2.2 Interactions and mixing with light fields 

Even though the MSSM fields are uncharged under the global symmetry, they may couple to 
the spontaneously broken sector through higher-order terms in the Kahler potential. We will 
particularly be interested in the coupling of the Goldstone multiplet with the Higgs doublets Hu,d- 
Explicit symmetry breaking terms can generate superpotential couplings between the MSSM and 
the Goldstone sector; these are discussed in Section |4j 

The Kahler potential interactions between the Higgses and the Goldstone superfield can be 
parameterized by expanding in 1//, 

K= ^{A + A^){c,H^Ha + ...+ h.c.) + ^{A + A^fic^H^H, + ... + h.c.) + 0{l/f) . (2.11) 

Note that the first term vanishes if there is a Z2 discrete symmetry A — > —A. The presence of 
such symmetry depends on the choice of UV completion. A mixing between the Higgs and the 
sGoldstone arises, for example, from the Kahler metric term 

K^f = d'K/{dH^dA^) = jC^H, + . . . ^ cos/3 + . . . . (2.12) 

The C2 terms can also give rise to mixing if the sGoldstone also gets a VEV of order (s) ~ /. 
After rotating the Higgs and sGoldtone fields, the coupling between h and the Goldstone multiplet 
appears in the effective Lagrangian as 



■eS 



1 i _ 

-{dado) + -xYdf,x 



b,^s + c,'^h + . ..]+..., (2.13) 
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where Ch is a function of the coefficients Ci_2, (ii,2 and the Higgs sector mixing angles. This couphng 
is suppressed in the large rus limit, Ch — > {mh/ms)'^. At this order in q^v-Em / P there are additional 
Higgs doublet couplings of the form 

' C2xYl'x{Hudf.H, + d,H^H,-h.c.) (2.14) 



4/ 

which give rise to additional interactions of the heavy Higgses with the Goldstone fermion, but 
do not involve the light higgs h. We neglect these couplings and the mixing of the heavy Higgses 
with the sGoldstone. 

Besides the scalar mixing, there is kinetic mixing between the Higgsino and Goldstone fermion 
of the form 



t-KM — ^ 



{x^cx^d^H^ - d^x^cr^Hu) (c^H, + ...) + h.c. + (H^ ^ H, 

^ ze„x^(T^9^if° + ^e,xV9^^° + h.c. (2.15) 

where eu,d ~ vew/ f- In the case where /i ^ m^, the Goldstone Fermion has a small Higgsino 
component roughly given by eu^m^/fi ~ vew^x/ ff^- 

The Kahler terms involving the other MSSM matter fields are typically more suppressed. 
Assuming minimal fiavor violation to control fiavor-changing neutral currents, these terms take 
the form 

K=hA + A^) ( ^QHuU+ ^QHdd+ ^LHde + h.c) . (2.16) 

The suppression scales Mu^d,i are not necessarily related to the global symmetry breaking scale /, 
and can be much larger depending on the UV completion of the theory. 



3 SUSY breaking 

We assume that soft SUSY breaking terms which simultaneously break the U(l) global symmetry 
are negligible. The remaining soft terms generate an explicit sGoldstone mass, but leave the 
Goldstone boson massless. The Goldstone fermion may only get a mass from the superpotential 
or from D-terms via mixing with gauginos. For simplicity we ignore the latter possibility so that 
the fermion mass matrix is the second derivative matrix of the superpotential, 

{mermion)ij = Wij. (3.1) 

While the superpotential terms are U(l) invariant and supersymmetric, a Goldstone fermion mass 
can still be induced if the vacuum is shifted from its supersymmetric value due to the presence of 
soft breaking terms. The U(l) invariance of the superpotential implies 



so that the Goldstone fermion x = Yli Qifii^i/f is indeed a zero mode of the fermion mass matrix 
when none of the U(l)-charged F-terms obtain a vev The mass of the Goldstone fermion 
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then depends on whether the U(l)-charged fields pick up F-terms in the presence of soft breaking 
terms [71 [8] : 

m^^qi{F,)/f. (3.3) 

If the superpotential has an unbroken i?-symmetry which is left unbroken by the soft terms, 
then X necessarily carries i?-charge —1 and a Majorana mass is prohibitecQ In particular, soft 
scalar masses always preserve i?-symmetry and hence cannot generate a Goldstone fermion mass 
in the i?-symmetric case. On the other hand, A-terms are holomorphic and generically break 
i?-symmetries. Thus A-terms are expected to always contribute to the Goldstone fermion mass, 
while soft scalar masses may or may not contribute. 

The effect of A-terms is equivalent to the mixing between the F-terms between U(l)-charged 
fields and the SUSY breaking spurion {X) = FO"^ + ■ ■ ■ . For concreteness, we consider gravity 
mediation with F/Mp\ ~ "^soft- It was recently emphasized in [TB] that FF^ type mixing terms 
are always expected and will contribute a mass of order 7/13/2 to the Goldstone fermion. Indeed, 
such mixing terms arise from higher dimensional Kahler terms of the form 

ir = ^Z(X,Xt)$J$i. (3.4) 

i 

Using the technique of analytic continuation into superspace [IT] , one may absorb Z into a redef- 
inition of the chiral superfields 

$ ^ = ^1/2 ( 1 + ^^^FeA $. (3.5) 



dX 

This canonically normalizes K and generates soft terms that include the A-terms 

dW 



AC 



soft 



9$ 



^-1/2 


d\nZ F\ 


dlnXMJ 





(3.6) 



These terms completely incorporate the mixing between F-terms of the form FF^^i. The Gold- 
stone fermion mass is determined by the induced FjS obtained by minimizing the scalar potential. 



V 



dW ^ 



+ A,|^0i + h.c. +m,2|0,|2. (3.7) 

0(pi 



To summarize this section, we find that A-terms will always contribute to the Goldstone 
fermion mass. Assuming that Ai,mi < fi for all the generic size of the induced F-terms is 
\Fi\ f» Aifi and, consequently, the induced Goldstone fermion mass is ~ AiQi. In many situations 
the A-terms can be suppressed relative to other soft breaking terms and it is thus reasonable 
to expect that the Goldstone fermion remains lighter than the other superpartners. Soft scalar 
masses may also contribute. If they do, their contribution to the Fj is expected to be of order 
Fj ~ mf so that the contribution to the Goldstone fermion mass is of order ~ mj/fi, which again 
can easily be suppressed. 



^We thank Y. Nomura for pointing out the role of i?-symnietry. 



7 



4 Super potential terms from explicit breaking 



The shift symmetry forbids any superpotential for the Goldstone chiral multiplet A. In order to 
generate a small Goldstone boson mass one must include terms which break the global symmetry. 
These can come from an anomaly in the global symmetry or through explicit breaking terms. 



4.1 Anomaly 

If the global symmetry is anomalous then the triangle diagram generates a aGG term which fits 
into a superpotential term 

W^anomaly = -yAiy'^VT'^ (4.1) 

where W"" = A'^ — ia^'^OG'^^^ + ... is the field strength chiral superfield for the gauge group G which 
has a U(1)G^ anomaly. In practice we will take G to be SU(3)coior or U(1)qed since we will be 
interested in the coupling to massless gauge fields. M4nomaiy generates non-derivative couplings in 
the effective Lagrangian: 

^anomaly D {aG^^G^ + ^X^^IJt^ (4.2) 

where G^^ = ^e^^p^G^^^. 

For the remainder of this document we assume that the global U(l) is anomalous. For example, 
the anomalous coupling Can is generated when the Goldstone boson a couples to fermions \E'j 
that transform in the fundamental of the gauged SU(A^) and carry a global charge qif,, 



57r 



/ f \ 



The result for a U(l) gauge group is similar and is obtained by including the different qi charges, 

c^J = '^V2N^j:2q^(y^), (4.4) 



an 



where Nc = 3 and the factor of 2 comes from the normalization of the generators in SU(A^), 
Tr [T'^T''] = 6''^/2. The simplest and most common case is when all masses are degenerate, 
mq,- = m^, and the yi are equal, yi = m^q,^ /{f\/2), so that 

a 

Can = ^q^N^ ■ (4.5) 

OTT 

Note that gauge coupling unification is preserved if the mediator fields \E'j are embedded in 
complete GUT multiplets. For example, one may consider A^^ x (5 © 5) representations of SU(5) 
which decompose intc[^ (3, 1)1/3 and (l,2)_i/2 under SUc(3) x SUl(2) x Uy(l). In this case the 
mediators are both colored and electrically charged; they thus allow the dominant decay to be 
a ^ gg with subdominant contributions from a — )■ 77 with branching ratio ~ 10^^. 



^The hypercharge normalization is fixed if there are no exotic electric charges. 
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4.2 Explicit breaking spurions 



Sources of explicit global symmetry breaking terms can be parametrized by spurion chiral super- 
fields Ra which carry charge a under the global symmetry and obtain a vev (Ra) = Aq,/, where 
\a ^ 1. This permits a superpotential term AW = pJ2a^-a^" where = exp{aA/ f). 
Unbroken supersymmetry requires that there are two sources of global symmetry breaking, R_a 
and R-/3, with opposite charges, a(3 < 0. This produces a sGoldstone boson vev and an effective 
superpotential with a common supersymmetry preserving masses = = m^. Explicit break- 
ing terms may also generate new interactions which are completely determined by the Goldstone 
boson mass, 

(4.6) 

The only model-dependent inputs are the charges a and (3 of the explicit breaking operators. 
After SUSY breaking, the sGoldstone boson and the Goldstone fermion masses are lifted while 
the Goldstone boson remains light, ^ rriy > ma- Up to integration by parts, the the on-shell 



trilinear axial coupling axi X is equivalent to an effective hi coupling in (2.7). Finally, the explicit 



breaking can also generate additional terms in the superpotential of the form 

W = cR^^G'^ (h^H^ + ^QH^u + ^QHd + ^LH^e + h.c.) . (4.7) 
\ Mu Md Ml / 

These lead to mixing with the Higgs and decays to SM fermions. 



5 Relic Abundance 

The Goldstone fermion x is a natural dark matter candidate if it is the LSP and produces the 
observed abundance [U HE] , 

f^D^/i^ = 0.112 ±0.0056, (5.1) 

where h is the Hubble constant. A key observation is that the effective interactions between the 
Goldstone fermion x and Goldstone boson a lead to an annihilation cross section xx ~^ '^'^ of the 
correct magnitude for a thermal relic with 0{1) couplings and mass at the SUSY breaking scale 

MsvSY, 



Note that an explicit factor of the temperature appears in (5.2) because parity forbids the s-wave 
channel. Thus the Goldstone fermion is an almost ideal WIMP candidate. Due to the slight 
thermal suppression, the coupling bi has to be slightly larger than 1. Otherwise, with the natural 
choices of parameters, the correct annihilation cross section is obtained. 
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5.1 Summary of model parameters 

Below we provide a summary of the Goldstone fermion model parameters and the values used in 
our parameter space scan: 



Parameter 


Description 


Scan Range 


/ 

ma 
bi 

Can 
Ch 

6= (/3-a)/2 


Global symmetry breaking scale 

Goldstone fermion mass 

Goldstone boson mass 

XXa coupling, (2.10) 

Anomaly coefficient, (4.2) 

Higgs coupling, (2.13 

Explicit breaking iaxi^X coupling, (4.6) 

Explicit breaking o?xX coupling, (4.6) 


500 GeV - 1.2 TeV 
50 - 150 GeV 
8 GeV - //lO 
[0,2] 
0.06 

[-1,1] 
3/2 

13/8 



These values represent a natural cross section of the full parameter space. 



5.2 Summary of annihilation channels 

In addition to 



(a) XX ~^ (^(^ ill the t-channel and w-channel via the self- interactions (2.7) 



a detailed analysis shows that there may be appreciable contributions from 



(b) XX — ^ '^'2 from explicit breaking terms (4.6) 



(c) XX — ^ 99 with a in the s-channel via the anomaly (4.2). 



In fact, these can overcome the p-wave suppression in the annihilation into 2a. Note that XX ~^ 99 
also gives an s-wave contribution which may contribute up to ~ 1/3 of the total annihilation cross 
section. Less significant are the decays into Higgs bosons. 



(d) XX — ^ with a in the s-channel via the Higgs coupling (2.13) when m/j -|- < 2m^; 



(e) XX hh via the coupling to two Higgs bosons (2.12) when nih < m^. 



Note that in some cases the Higgs boson may be lighter than the 115 GeV because of non-standard 
Higgs decays (see Section[7]for the relevant collider phenomenology). Other annihilations involving 
a virtual gluino, the sGoldstone, or the Higgs boson are typically suppressed by large masses or 
small Yukawa couplings. A detailed calculation of each contribution is presented in Appendix |B} 
The model generates the correct abundance for Goldstone fermion masses between 50 — 150 GeV 
and Goldstone boson masses between 10% — 100% of for couplings bi ~ 0{1). Fig. [Tjshows the 
contours for different values of rria/m^ subject to the correct relic abundance in the {m^, bi) plane. 
It may further be possible to open up a different region of parameter space with lighter Goldstone 
boson and fermion scales through a Sommerfeld enhancement due to an attractive force between 
the exchange of multiple low-energy Goldstone bosons [H]. We briefly discuss this possibility in 
Appendix [O 
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Amines havenA-=0.11 




60 80 100 120 140 

/H^[GeVl 



Figure 1: Contours for different values of the Goldstone boson mass: m^/m^ = 0.1 (black dotted), 
0.5 (blue dashed), and 0.7 (red solid) for fixed relic density Qh^ = 0.11 in the (m^, bi) plane. Gray 
lines include the subleading contributions from annihilations into Higgs bosons, XX ~^ and 
XX ~^ hh. The kink at 60 GeV comes from threshold effects. We set / = 700 GeV, a = —4, 
/3 = 1, m/, = 116 GeV, Can = 0.06. 



6 Direct and Indirect Detection 

We have seen above that this model can easily produce the correct dark matter abundance. Next 
we estimate the generic size of direct and indirect detection bounds. 

6.1 Dark Matter Effective Operators for Direct detection 

In order to evaluate the cross section of Goldstone fermion scattering off nuclei in direct detection 
experiments, one must evaluate the nucleon matrix elements. Usually one parametrizes the light 
quark mass content fi^^ of the nucleons, 

mi{N\qiqi\N) = f^^^ruN i = u,d,s, N = p,n. (6.1) 

where m^v is the nucleon mass. The largest contribution comes from the strange quark [20], 
but with sizeable uncertainties |2T]. We assume the default value in the micrOMEGAs code, 
fifd ^ ~ fs^^ ~ 0-26 [22] • For heavy quarks, the contribution fjf^'' is induced via gluon 
exchange and can be calculated by means of the conformal anomaly |23j . 

m,(iV|g,g,|iV)^/fW = ^m^M - /f M , h = cAt. (6.2) 

\ i=u,d,s / 
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6.1.1 Coupling to quarks via Higgs exchange 



In Section |2.2| we showed that after electroweak symmetry breaking, it is natural to expect a 
non-vanishing coupling between x aiid the lightest neutral Higgs boson /i, 

Ch = Ch^^{xiYdf,x)h, (6.3) 

where the size of the coupling Ch depends on the specific realization. The Higgs coupling to 
nucleons is set by the Yukawa couplings and — in the presence of more Higgses — the mixing angles, 
CqTUq/ {\/2vEw)hqq- Integrating out the Higgs generates an effective four-Fermi interaction, 

= G^.NNxx G,, = c,0= (^) , (6.4) 



where we used the equations of motion for x and the quark content of the nucleons (6.2) to write 



q=u,d,s \ q=u,d,s / \q'=c,b,t / 

The resulting scattering cross section per nucleon at zero momentum transfer i^ 

afr=^[G^pZ + G,„{A-Z)]\ (6.6) 

where /i-^ = {m~^ + m^^)^^ is the reduced mass. The typical value for crf^^^^ is just below the 
XENONIOO direct detection bound 

. S4 X :o-« ..^ {^)' (mf^y (^J , 

Note the {m^v/PY suppression present in this cross section (due to the Goldstone nature of x) 
relative to that of a generic Higgs exchange. For example, Higgs-mediated neutralino decay in 
the MSSM with coupling C ~ cg/2xxh needs a very small coupling c to avoid the XENONIOO 
bounds: 

27r mlvl^ 

Thus, Goldstone fermion dark matter offers a natural suppression of the direct detection cross 
section while retaining the correct WIMP annihilation cross section and abundance. Fig. [2] plots 
typical values of the direct detection cross section for parameters with correct relic abundance. 

^ Since most direct detection events occur at low recoil energy, it is standard to parameterize the cross section in 
terms of a zero momentum transfer part and a form factor which encodes the momentum and target dependence. 
See, for example, |24) . 
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Figure 2: Black line: XENONIOO bound. Left: scan over parameter space with 500 < / < 
700 GeV (blue), 700 </ < 800 GeV (violet), 800 </ < 900 GeV (green), 900 </ < 1000 GeV 
(yellow). We scan < 6i < 2, -1 < Ch,hh < 1, 8 GeV < m„ < ^, 50 GeV < < 200 GeV. 
Right: Blue points have 0.5 < nia/m^ whereas red points have ma/rriy < 0.5. 



6.1.2 Coupling to gluons 



Integrating out the massive gaugino in (4.2) generates two dimension-7 operators that couple x 
to gauge bosons, 

= - (^) [XX] Gl,Gl, = -^ [xi'x] Gl.Gl, , (6.9) 

where Mx is the gaugino mass. In the limit of zero momentum transfer, only contributes 
to direct detection since GG is a total derivative. We therefore neglect C^^j hereafter. The 



df} can be mapped to a standard four-Fermi operator 



{N\GG\N) nucleon matrix element can be extracted from the conformal anomaly (6.2) so that 

-eff 




4^ 4ff!^ = G^XXNN , = I 1 - > : fr' I • (6.10) 

The corresponding cross section per nucleon at zero momentum transfer is 

2 /^7rin r-i Tr\ 4 



_93 
'SI 



----^mm(^)(f)'(riv)^ (a.n, 

where Can = agq^iNxi, / {St[) has been used. This value is much smaller than both the recent 
upper bound by the XENONIOO collaboration [3j and the expected reach at the LHC, = 
few X 10-^6 cm2 



6.2 Indirect detection 

Many experiments are searching for indirect signals of annihilation of dark matter in dense en- 
vironments such as the galactic center or the solar core. The rate of such events is set by the 
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present-day thermally averaged annihilation cross section. Note, however, that the dominant an- 
nihilation channels at freeze-out are p-wave and hence are strongly velocity suppressed in the 
current era. Thus, the relevant annihilation channels for indirect detection are s-wave and were 
sub-dominant at freeze-out. These cross sections are relatively small and astrophysical observa- 
tions do not impose severe constraints. 



6.2.1 Fermi-LAT: lines, isotropic diffuse 7-rays, and dwarf galaxies 

Dark annihilation in the galactic halo may produce photons either directly (e.g. XX ~^ ll) 
through secondaries (bremsstrahlung off charged products or decays of neutral pions). The Fermi 
experiment has searched for excesses in the gamma ray spectrum both in the form of lines arising 
from prompt annihilation to photons and in contributions to the diffuse spectrum from secondary 
products of annihilation. 

Fermi currently searches for 7-ray lines from 30 - 200 GeV PH], with upcoming bounds that 
are an order of magnitude stronger in the 7-30 GeV region [27]. The lack of a bump in the 
Fermi data implies an upper bound on (crf)^^ between (0.2 — 2.5) ■ 10~^'^cm'^/s when using the 
Einasto dark matter halo profile which predicts the largest photon flux among those examined in 
the Fermi analysis. 

In the Goldstone fermion model, prompt annihilation to photons occurs through an anomaly 
vertex similar to the one which mediates a — ?■ gg. This rate depends on the U(l) x U(l)|f^ anomaly 
coefficient which is determined by the choice of electric charges for fields carrying global charge. 
The cross section for annihilation into gluons is given in ( B.l[ ). The analagous expression for 



annihilations into photons is given by replacing alN^ — )■ a^-^ {'^Yliii.Qku)'^T ■ the case where 
the ^ are taken to be in the (anti-)fundamental of an SU(5) unified group, we find (cv)^^ ~ 
2 ■ 10^^ {av)gg. Even with the most extreme choices of the model's free parameters, this rate 
remains more than an order of magnitude smaller than the Fermi bounds. 

Fermi has also measured the isotropic diffuse 7-ray spectrum in the range 20 — 100 GeV 
|28j . This bounds the annihilation of dark matter into charged particles and neutral pions. For 
example, for a 400 GeV dark matter particle which annihilates into a bb pair, Fermi sets a bound 
on {crv) f^i which is roughly an order of magnitude above the cross section required to reproduce the 
right relic abundance. The Goldstone fermion model generates diffuse photons primarily through 
annihilation to gluons produced in the s-wave annihilation channel xX ~^ 99- However, this cross 
section is at least an order of magnitude smaller than the Fermi bound and hence the Fermi diffuse 
7-ray data do not constrain this model. 

Preliminary results from a Fermi analysis of 10 dark-matter-rich dwarf spheroidal galaxies also 
place limits on photo-production from dark matter annihilation [22] ■ For low-mass (< 60 GeV) 
dark matter annihilating into bb pairs, constraints on the annihilation rate extend slightly below 
the thermal relic rate of 3 ■ 10~^^ cm^/s, with the strongest constraint of ~ 1 ■ 10~^^ cm^/s at 
= 10 GeV. In this mass window and for reasonable parameter choices, the Goldstone fermion 
annihilation cross section is always at least a factor of 3 lower than these limits. 

Other constraints, such as those that come from 7-rays originating in clusters of galaxies, 
typically set weaker bounds 
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6.2.2 PAMELA: the antiproton flux 



Antiproton flux Antiproton flux 




Kinetic energy[GeV] Kinetic enei"gy[GeV] 



Figure 3: Antiproton flux at the Earth for / = 700 GeV, = 2, 6 = 3/2, N^;, = 5 at fixed density 
Qh"^ ~ 0.1. Red (blue) lines represent the propagation parameters MAX (MIN) used in [33] with 
the Einasto DM halo profile. The dots represent the PAMELA data [3T]. Left: nia/m^ = 0.5 at 
= 50 GeV and 6i = 3 (solid); = 100 GeV and bi = 1.5 (dashed); = 150 GeV and 
bi = 1 (dotted). Right: nia/m^ = 0.8 at = 50 GeV and bi = 2.5 (solid); = 100 GeV and 
bi = 1.2 (dashed); = 150 GeV and bi = 0.5(dotted). 



PAMELA has recently published data on the absolute cosmic ray antiproton flux from 60 MeV 
- 180 GeV ^31J. This places constraints on dark matter models with a substantial annihilation 
rate to hadrons. For a 100 GeV WIMP, the annihilation cross section to Zs, Ws, and b quarks 
has an upper bound comparable to the rate required for the observed relic abundance, (cf )reiic ~ 
3 ■ 10~^^cm^/s [32]. For Goldstone fermion dark matter, the dominant annihilation channel in 
the galactic halo, XX ~^ 99: is s-wave. This has a typical cross section of (av) ~ 10~^^ cm^/s 
and can be pushed up as high as 10~^^ cm^/s. Using recent numerical recipes [33], one may 
estimate the anti-proton flux as a function of the thermally averaged annihilation cross section 
and the Goldstone fermion mass. This is depicted in Fig. [3] for different model parameters that 
yield the correct relic abundance. The anti-proton flux varies considerably as a function of the 
galactic propagation parameters and the halo profile. The solid, dashed, and dotted curves each 
correspond to different underlying Goldstone fermion model parameters. Choosing different halo 
profiles and propagation parameters leads to a spread in the predicted p flux such that for each 
choice of Goldstone fermion model parameters, the actual flux from dark matter annihilation is 
expected to lie between the two solid, dashed, or dotted curves respectively. 

For each choice of model parameters, there is a sizeable region where the predicted flux from 
dark matter annihilation lies well below the measured anti-proton flux. Thus the PAMELA data 
do not place signiflcant constraints on the Goldstone fermion dark matter model. 
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7 Collider Phenomenology 



In addition to (in) direct detection, Goldstone fermion models lend themselves to novel collider 
signatures coming from the Goldstone supermultiplet. As discussed in Section [3| the sGoldstone s 
is typically heavy with small couplings to the SM sector so we may neglect its collider signatures. 

7.1 Collider signals of dark matter 

The most direct way of testing the dark matter annihilation mechanism is through dark matter 



pair production coming, for example, from the XX^G operator in (6.9). One signature of this 
process at colliders is a monojet coming from hard initial state QCD radiation. For a range of 
masses up to the TeV scale, the LHC will set the most stringent bound on this operator with a 
sensitivity of cxg^ ~ 10~^^ — 10~^^cm^ for a 5a discovery with 100 fb~^ ^25]. The effective scale 



that suppresses the dimension-7 operator (6.9) is roughly an order of magnitude larger than the 
LHC reach. However, the process gg ^ a* ^ XX via an off-shell Goldstone boson gives a larger 
contribution and results in a naive effective scale M* ~ {mxP / '^anY^^ ~ 1 TeV. This is in the 
ballpark of the LHC 5a reach given in |25] . 

Goldstone fermion dark matter can also be produced from the cascade decay of heavier i?-parity 
odd particles, such as gluinos or squarks. Due to the small coupling between the MSSM and the 
Goldstone sectors, the cascade decays will all go through the lightest ordinary supersymmetric 
particle (LOSP). The decay of LOSP to the Goldstone fermion is determined by the operators 
connecting the two sectors. In the current setup, there are two types of interactions: 



the anomaly induced coupling xGX, as in (4.2), and 



the kinetic mixing discussed in Section 2.2 



The details of the decay modes depends on the nature of the LSP. For example, a bino-like LOSP 
will decay to the LSP via the anomaly, B — )■ x + l/Z. A Higgsino-like LOSP would decay instead 
to the LSP because of the kinetic mixing, /i — x + /i, and — ?■ x + ~^ X + 2j- In the latter 
decay mode, the reconstruction of the Goldstone boson resonance in the jet final state is difficult 
if a is below 100 GeV, but it may be possible instead in the diphoton decays of a with sufficient 
luminosity. These channels yield prompt decays even though they may be suppressed by loops or 
small mixing angles. For example, the natural width for a pure bino LOSP is around 10^^ GeV. 

Finally, the presence of exotic heavy fermions also has interesting implications at colliders. 
These fermions may be considered to be "fourth generation" quarks which, if they are sufficiently 
light, can be probed at the early stages of the LHC (see the discussion in [T3] for an example). 

7.2 Non-standard Higgs boson decays 

The largest natural coupling of the Goldstone boson and fermion to the SM is through Higgs 



boson via the kinetic terms, (2.13). This coupling allows the Higgs to decay into 2a or 2% if 
kinematically allowed. Typical branching ratios are plotted in Fig. |4j 

The Higgs boson decay h ^ 2a gives rise to four light, unflavored jets coming from a 2g. 
This decay mode is therefore easily 'buried' under the QCD background. Such non-standard 
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Figure 4: Plots of Higgs boson branching ratios for various parameters. 



Higgs boson decays have recently been investigated in SUSY models where the Higgs boson itself 
is also a pseudo-Goldstone boson emerging from the spontaneous breaking of a global symmetrj]^ 
In particular, the spontaneous breaking of SU(3) — )• SU(2) gives rise both to a light Goldstone 
multiplet A and a light Higgs multiplet [Ej. The resulting coupling Ch ~ \/2 is set by the kinetic 
mixing between the two multiplets which, in turn, is fixed by the scale / of the global symmetry 
breaking. A more recent example of a 'buried Higgs' in SUSY has been discussed in the context 
of a spontaneously broken U(l) symmetry where Ch depends on couplings in the superpotential 
because the Higgs is no longer a pseudo-Goldstone boson [15] . 

Even though these non-standard Goldstone fermion decay modes can dominate, the branching 
ratio to SM particles is still larger than ~ 20% at low Higgs masses and therefore the LEP bound 
on the Higgs mass cannot be lowered below ~ 110 GeV. Furthermore, while the discovery of a 
completely buried Higgs is challenging at the LHC [36], this 'partially buried' Higgs would be 
discovered in SM channels with a missing piece in the total width. The invisible Higgs boson 
decays (xs leave the detectors) can be probed at the LHC through the missing energy signal [37] . 
Both the buried and invisible decay modes may have sizeable branching ratios, and the observation 

^For early attempts of this idea in SUSY see [34]. More recently, SUSY and little Higgs models motivated by 
the little hierarchy problem have been proposed |35) . 
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of both channels would give strong evidence for this scenario. 

8 Conclusions 

Acceptable dark matter scenarios within the MSSM must become increasingly contrived as the 
sensitivity of direct detection experiments increases. In order to remain consistent with recent 
XENONIOO results, neutralino WIMP models must typically invoke accidental mass relations to 
boost the annihilation cross-section through co- annihilations or strategically placed resonances. 

Inspired by this tension, we have explored a general supersymmetric framework compatible 
with GUT unification in which the LSP is the fermionic component x of a Goldstone supermultiplct 
associated with a U(l) global symmetry that is spontaneously broken at the TcV scale. Because 
the Goldstone fermion's couplings to the Standard Model are suppressed by ~ ^'Ew^x//^ 
additional loop factors in some cases) , these models are able to avoid direct detection constraints 
from XENONIOO and indirect detection constraints from Fermi and PAMELA. 

The annihilation cross section of a weak-scale Goldstone fermion at freeze out is on the order 
of 1 pb, with dominant contributions coming from p-wave annihilation into Goldstone bosons. 
Typically subdominant s-wave annihilations into gluons arise through anomalies of the new global 
symmetry. The observed dark matter relic density is obtained with natural values for the model 
parameters. 

This class of models also offers novel and distinctive signatures at colliders. Goldstone fermions 
can be produced at the LHG in pairs through the anomalous coupling to gluons, leading to monojet 
signals when there is additional hard QCD radiation from the initial state. Additionally, SUSY 
cascades are modified by decays of the NLSP to the Goldstone fermion. Examples include the 
bino decay to a photon and the Goldstone fermion, and the higgsino decay to the Goldstone 
fermion and the Goldstone boson. The Goldstone multiplet also modifies the phenomenology of 
the Higgs sector. Interactions with the Goldstone boson allow cascade decays of the Higgs to four 
jets, h ^ 2a ^ 4:j, analogous to models where the Higgs decays are 'buried' under the QCD 
background. If kinematically allowed, the Higgs may also have a sizeable fraction of 'invisible' 
decays, h XX- 
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A Explicit models 



We present explicit models to demonstrate how one may generate different values of the coupling 
bi, defined in (2.10). In their simplest form, both examples have an unbroken i?-symmetry which 
implies that only the A-terms generate a mass for the Goldstone fermion. It is straightforward 
to modify these examples to explicitly break the i?-symmetry without modifying the structure of 
these theories. 



A.l The simplest example 

We consider a simple variation of the model considered in [TB] with the superpotential W = 
yS{NN - 12^). This gives 

K = fle^^^^'yf + /le-(^+^^)// f = fN + fN, (A.1) 

so that the tree-level range for bi is 



Pn - P 
In ~ In 



-l<bi = '-^ ^ < 1 . (A.2) 



A. 2 An example with |6i| > 1 

A perturbative model that may give > 1 is the following: 

W = \XYZ - ^i^Z +^Y'^N -^NN , (A.3) 

where the charges are qz = 0, = —Qn = ~2gy = 2qx and all couplings and masses are non-zero. 
The resulting supersymmetric minimum 

fxfY = fiVX, fz = fN = 0, fN = ~X^ (A.4) 

2/i 

gives vanishing F-terms while the Goldstone chiral multiplet is 



f ~ f 

The corresponding bi at tree-level is given by 



which goes to 6i — )■ 2 when ffq ^ fx,Y- 
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(a) XX 99 



(b) XX aa 



(c) XX aa 




X h 




X a 
(d) XX aa (e) XX ^ ha (f ) xx ^ hh 

Figure 5: Goldstone fermion annihilation channels. 

B Annihilation cross section 

Diagrams for the dominant annihilation channels are presented in Fig. [5j 
B.l XX 99 



The annihilation cross section to gluons (see Fig. 5a) is controlled by the anomalous coupling (4.2) 
where Can = asq^!iNq,/{8TT) and the vertex 6i/(2A/2/)x7^7^X(9^a. Away from resonance one finds 



av 



where v is the relative velocity in the center of mass frame, A'^c = 8 is the number of colors in the 



final state, and 5 = —{a + /3)/2 is the contribution from the explicit breaking vertex (4.6). Note 
that this process gives a non-vanishing s— wave annihilation component. 



B.2 XX aa 



Annihilation into Goldstone bosons proceeds through t- and u- channel diagrams (see Figs. 5b -5c) 
as well as a contact interaction coming from explicit breaking (see Fig. 5d). 
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B.2.1 t- and w-channel 

These diagrams give a p-wave contribution to the cross section, av = a + hv"^ + . . ., 

a =0 z = nia/m^ (B.2) 

2 

^ = na fA( 2 o^4 [^i(^i + 4^'^)(3^' - 16^' + 48^^ - 64^2 + 32) (B.3) 



where 5 = — (a + /3)/2 is the contribution from the exphcit breaking vertex (4.6). 



B.2. 2 Explicit breaking vertex 

The quartic contribution to the annihilation cross section is also p-wave is 



0"^ = pSrr^a -\ Va = \l (B.4) 



where p is given in terms of the charges of the explicit breaking operators (4.6), p = a + a(3 + f3 



B.2.3 Interference 

The contact interaction interferes with the t- and u-channel diagrams. Summing the amplitudes 
and then squaring gives, 

^ = ?^(26?+86i^5+^p) + — ^ (3p2 + 32h6p + 128bl6^ - 166^) +0{z^) , a = (B.5) 
y67r/^ i5367r/^ 

where av = a + bv^ + . . . and z = ma/my.. Note that for all plots in this document we use the full 
expression for b that is valid for all z < 1. 

B.3 Subleading processes 



The annihilations to a single Higgs (Fig. 5e) and to two Higgses (Fig. 5f) are subdominant. 



B.3.1 XX ^ ct* — ^ 

This channel is available when 2m-^ > rria + ruh. Naively, it should be less important because 
the cross section has an extra suppression by (wew//)^- On the other hand, this is an s-wave 
contribution and therefore the effect is not completely negligible compared to the xx aa p- 
wave process. The cross section is given by 

Va [bim^ + 5mafclvl^ fml-ml + s 



a 
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B.3.2 XX ^ hh 



This channel is allowed when > mh, up to thermal contributions. Because Higgs can be 
buried under QCD, it is possible to have rrih ~ 90 GeV. This process is generated from the 
contact interaction term c/i/j(x^7^<9^x)^^/(2/^), which follows from the C2,d2 coefficients in the 
Kahler potential. 



,2^2 /„ 4^2 



87r / 

This is a p-wave process. 



av = ^'^[- ^) (B.7) 



C Sommerfeld enhancement from Goldstone boson ex- 
change 

Thus far we have calculated the relic density assuming no enhancement due to long-range forces. 
Here we briefly present the non-relativistic potential between the Goldstone fermions and argue 
that there could be regions of parameter space with a sizeable Sommerfeld enhancement in the 
annihilation cross section due to an attractive force between the Goldstone fermions due to the 
exchange of multiple low-energy Goldstone bosons [19] , as depicted in Fig. [6] It is thereby possible 
to lower the Goldstone boson and fermion mass scales. We emphasize that this enhancement is 
not necessary to obtain the correct abundance and sufficiently low direct detection cross sections, 
but it may open up a different region of the parameter space where the Goldstone fermion mass 
in the 10 — 50 GeV range. 



Figure 6: Exchange of multiple soft Goldstone bosons can lead to an attractive force enhancing 
the annihilation cross section for the Goldstone fermions. 



In the non-relativistic limit, the X1X2 — ^ XvX2> scattering amplitude gives rise to a spin-spin 
interaction. The low-energy potential can be written in terms of a traceless tensor and a central 
piece: 

V{r) = Vrir) (s Si ■ f S2 ■ f - Si ■ S2) + Vb(r)5i • ^2 , (C.l) 
where the coefficients are 

V.,M = JL(i + !!^+l!!!l)e— Vc(r) =JLl!!ie— (0.2) 

Svr/^ \r-^ 6 r J Snj^S r 

Note that the leading term for distances r < is contained in the tensor potential. For total 
spin S = 0, the tensor potential averages out to zero, whereas the central part gives an attractive 
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interaction which is independent of the orbital angular momentum 



{S = 0,i\V{r)\S = 0,i) = -^Vc{r). 



(C.3) 



This contribution vanishes in the limit — )■ in agreement with |3H]- For S = 1, i = 1 the central 
potental is repulsive whereas the tensor is attractive. The net effect is an attractive potentialj 




The magnitude of this Sommerfeld enhancement was calculated in detail for s-wave annihilation 
processes in [12], where it was found to take values as large as 1000 and as small as 0.1. For 
the current model one would only need a factor of few to lower the Goldstone boson and fermion 
masses to the 10 GeV range. Since most of the leading annihilation channels relevant to this class 
of models are p-wave, the results of [19] are not directly applicable. A dedicated calculation is left 
for future work. 
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